Quantum Simulation of Spin Models on an Arbitrary Lattice with Trapped Ions 
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A collection of trapped atomic ions represents one of the most attractive platforms for the quantum 
simulation of interacting spin networks and quantum magnetism. Spin-dependent optical dipole 
forces applied to an ion crystal create long-range effective spin-spin interactions and allow the 
simulation of spin Hamiltonians that possess nontrivial phases and dynamics. Here we show how 
appropriate design of laser fields can provide for arbitrary multidimensional spin-spin interaction 
graphs even for the case of a linear spatial array of ions. This scheme uses currently existing trap 
technology and is scalable to levels where classical methods of simulation are intractable. 



PACS numbers: 

Quantum systems are notoriously difficult to model ef- 
ficiently using classical computers, owing to the exponen- 
tial complexity in describing a general quantum state as 
the system grows in size. In the 1980s, Richard Feynman 
proposed to circumvent this problem by employing a con- 
trol quantum system with tailored interactions and logic 
gates between quantum bits (qubits) in order to simulate 
the quantum system under investigation [1, 2]. 

Quantum spin models such as the Ising model have 
become a proving ground for Feynman's proposal, with 
systems of qubits behaving as effective spins and inter- 
actions engineered with external electromagnetic fields. 
The quantum Ising model is the simplest spin Hamil- 
tonian that exhibits nontrivial aspects of quantum mag- 
netism such as spin frustration, phase transitions [3] , and 
poorly understood spin glass and spin liquid phases [4, 5]. 
Indeed, solving for the ground state configuration of spins 
subject to a general fully-connected Ising interaction is 
known to be an NP-complete problem [6] . Here we show 
how a fully-connected Ising or more general Heisenberg 
spin model with arbitrary couplings across the spin net- 
work can be generated in a scalable system of trapped 
atomic ions, even for a one-dimensional chain in space. 
This may allow quantum simulations with hundreds of 
spins, where the physics cannot generally be predicted 
otherwise. 

Isolated atoms are an ideal control quantum system for 
quantum simulations, as they exhibit very long coherence 
times, can be measured with high efficiency, and can ob- 
viously be replicated with nearly identical characteristics. 
There has been great progress in the use of atoms trapped 
in optical lattices for studies of transport phenomena 
such as those described by Hubbard models [7] and the 
fractional quantum Hall effect [8] . Following original pro- 
posals to exploit atomic ions for the simulation of spin 
models [9-11], there has been steady progress in simu- 
lating Ising models with an effective transverse magnetic 



field, showing the behavior of quantum phase transitions 
in the limit of less than ten spins [12-18]. Recently, 
the Ising model with an axial field was demonstrated in 
neutral atoms in optical lattices with nearest-neighbor 
couplings [19]. All of these studies involve global spin- 
dependent optical dipole forces to generate trivial forms 
of the spin couplings, such as a uniform ferromagnet. Be- 
low we show how to tailor optical forces to generate arbi- 
trary fully-connected networks of N spins that uniquely 
specify each of the N(N—1) /2 pairwise interactions. The 
scheme is independent of the spatial geometry of the ion 
crystal and is compatible with one-dimensional arrays of 
trapped ions used in current experiments. 

We start with the arbitrary fully-connected Ising 
Hamiltonian on N spins, 

// y^.i,.,~; a/ . (1) 

where the matrix Jij is the strength of the spin-spin 
coupling between atoms i and j. The Pauli spin operator 
ai*' refers to the effective spin- 1/2 system within each 
atom, here represented for example by a pair of hyperfine 
ground states separated by frequency uj s [20]. 

The spins are coherently manipulated through a pair 
of counter-propagating laser beams that drive stimulated 
Raman transitions between the spin states while also cou- 
pling off-resonantly to the collective motion of the atomic 
chain [21, 22]. The atoms are arranged in a linear array, 
as is typical in a linear radiofrequency ion trap, although 
this scheme can also apply to other spatial geometries 
[11, 23, 24]. When the difference frequency between the 
Raman fields is bichromatic, with two spectral compo- 
nents tuned symmetrically at uj s ± fi with /i <C lj 8 , the 
effective spin-spin interaction of Eq. 1 emerges, medi- 
ated by the Coulomb-coupled motion of the atomic ions 
crystal [25, 26]. The axis of the spin-spin interaction on 
the Bloch sphere can be controlled by the relative optical 
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FIG. 1: Spectrum of transverse mode frequencies w m for 
AT = 10 ions in an anisotropic harmonic linear trap (solid 
lines), with the highest mode frequency uj\ corresponding to 
center-of-mass motion. Raman beatnote detunings /x m from 
the qubit frequency uj s are denoted by the TV = 10 dashed 
lines, with each spectral feature near a given motional side- 
band. The height of the dashed lines represents the intensity 
of each beatnote for ion i. In general each ion will be illumi- 
nated with a different set of intensities. 



phase of the two Raman fields [27, 28]. We assume the 
Raman lasers have wave vector difference 8k along the 
principal X— axis of transverse motion of the ion crys- 
tal [29] so that each X— transverse normal mode m with 
frequency uj m contributes to each spin-spin coupling [13], 
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Here, the Lamb-Dicke parameter r\i, m = 
bi, m fiky' n/ '(2Mw m ) sets the scale for the coupling 
between spin i and mode m, b i m is the normal mode 
transformation matrix, M is the mass of a single atom, 
and the single spin Rabi frequency fij is proportional to 
the optical intensity at atom i. The above expression 
is valid when the transverse modes are all laser-cooled 
to the Lamb-Dicke limit [21, 30], which is routine in 
typical linear traps [13]. We also assume the symmetric 
detuning /i is set sufficiently far from all motional 
sidebands (\u) m — fi\ >> r?i, m fii), so that the phonon 
states can be adiabatically eliminated, leaving the pure 
spin-spin coupling above [13, 25, 26]. 

The above expression has N + 1 control parameters in 
the set of Rabi frequencies {f2i} and the global beatnote 
detuning fj,. In order to generate an arbitrary Ising cou- 
pling matrix Jjj however, it is necessary to have at least 
N(N—1) /2 independent controls. Additional control pa- 
rameters can be introduced by adding multiple spectral 
beatnote detunings to the Raman beams, one near each 
motional mode (see Fig. 1), with a unique pattern of 
spectral components on each ion. There are several ways 
to achieve this, all involving some form of individual ion 
addressing. For simplicity, we retain the same set of N 
Raman beatnote detunings {/x m } on each ion and allow 
the spectral amplitude pattern to vary between ions, all 
characterized by the N x N Rabi frequency matrix f2i : „ 
of spectral component n at ion i. Note that the rela- 
tive signs of the Rabi frequency matrix elements can be 
controlled by adjusting the phase of each spectral com- 
ponent. This individual spectral amplitude addressing 



provides TV 2 control parameters, and the general Ising 
coupling matrix becomes 
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where Fij n characterizes the response of Ising coupling 
Ji j to spectral component n. An exact derivation of the 
effective Hamiltonian given a spectrum of spin-dependent 
forces gives rise to new off-resonant cross terms, which 
can be shown to be negligible in the rotating wave ap- 
proximation, as long as the detunings are chosen so that 
their sums and differences do not directly encroach any 
sideband features in the motional spectrum of the crystal 
[29]. 

We tune each beat note frequency near a unique nor- 
mal mode so that -Fi.^n has independent contributions for 
each n. Given a desired Ising coupling matrix Jjj, we use 
standard constrained nonlinear optimization to find the 
corresponding Rabi frequency matrix £li tU , while mini- 
mizing the total beam intensity. The deviation between 
the desired and the attained coupling was less than typ- 
ical round-off errors. 

We now present two example solutions for fij >n that 
produce interesting interaction graph topologies. First 
we calculate a Rabi frequency matrix that results in a 2D 
square lattice of nearest-neighbor antiferromagnetic cou- 
plings with N = 25 ions (5x5 grid with periodic bound- 
ary conditions), shown in Figs. 2a-b. Next we produce 
a 2D Kagome lattice of antiferromagnetic interactions, 
a geometry that can support high levels of geometrical 
frustration [31], shown in Figs. 2d-e. In both cases we 
assume the center-of-mass (COM) mode to be lui/2tt = 5 
MHz , and a fixed total optical intensity corresponding to 
Yli n \^i,n\ = 1 MHz . The beatnote frequencies fj, m are 
each tuned blue of the mode m sideband by a fraction f s 
of the spacing u>i — u>2 between the most closely-spaced 
modes (the COM and "tilt" modes, see Fig. 1), which 
itself scales as log/V/TV 2 . In these examples, the sparse 
nearest-neighbor nature of the interaction graphs require 
that most of the Ising interactions vanish, indicating a 
high level of coherent control over all of the Ising cou- 
plings. 

In order to generate a unique spectrum of Raman 
beams for each of N ions, some type of individual ad- 
dressing is necessary. For simplicity, we assume one of the 
two Raman beams is uniform and monochromatic, and 
the high frequency beatnote near the qubit frequency cu s 
or other global offset frequencies can be set by tuning this 
monochromatic beam. We focus attention on providing 
the requisite frequencies of the second beam, spread over 
a range given by the bandwidth of the transverse mo- 
tional mode frequencies of the ion chain, typically in the 
range 1 — 5 MHz . We suggest three possible methods for 
providing spatial dependent frequency modulation to one 
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FIG. 2: (a) Calculated Rabi frequency matrix Oi,„ to gener- 
ate 2D square lattice shown in (b), using the linear chain of 
TV = 25 ions shown in (c). The ion index refers to the order 
in the linear chain. The attained Jij nearest-neighbor is 27.6 
Hz for f s = 0.1. (d) Calculated Rabi frequency matrix fij jn to 
generate 2D Kagome lattice shown in (e) using a linear chain 
of TV = 36 ions. The attained Jij nearest-neighbor is 93.4 Hz 
for f s = 0.03. In both cases the total optical intensity corre- 
sponds to a Rabi frequency of 1 MHz if focused on a single 
ion, the nearest-neighbor couplings are antiferromagnetic and 
we impose periodic boundary conditions. 



of the Raman beams. The first method (Fig. 3) splits a 
single beam with a linear chain of TV individual optical 
modulators (e.g., acoustooptic or electrooptic devices), 
driven by TV independent arbitrary waveform generators. 
The second method splits a single monochromatic beam 
into a, NxN square grid and directs them onto a 2D array 
of TV 2 micromirrors [32] that are each individually phase 
modulated at a single frequency (and phase) [33] and fi- 
nally focused on the ion chain. The third method again 
splits the beam into an TV x TV grid of beams, this time 
with the vertical direction split by a single acoustooptic 
modulator, correlating beam position to frequency. This 
beam is then directed into a spatial light modulator that 
acts to mask (or phase shift) each of the TV x TV beams 
independently, and again focused onto the ion chain. In 
these implementations, it may be desirable to work with 
a uniformly spaced array of ions in the linear trap, so 
that the modulating elements arc also uniformly spaced. 
This can be accomplished by using a quartic or higher 
order linear trap [34, 35]. 

As the number of spins TV grows, the optical mod- 
ulation scheme becomes more complex, with either N 
or TV 2 elements required. We now estimate how the 
Ising couplings are expected to scale with the number 




FIG. 3: Schematic for individual spectral addressing a linear 
chain of N ions. A laser beam is split into a linear array of 
spots that each traverse TV independent acousto- or electro- 
optical modulators, driven by an TV independent arbitrary 
waveform generators (AWG). Alternatively, as discussed in 
the text, the beam can be broken into an array of N 2 beams 
that strike a TV x TV array of micromirrors each independently 
modulated, or a spatial light modulator. 



of spins along with errors due to experimental fluctua- 
tions, phonon creation and spontaneous emission scat- 
tering, assuming a fixed transverse mode bandwidth. 
The probability of phonon creation scales as p p h — 

Yli m ( ju'n-itm ^) ' ^ e off-resonant optical dipole forces 
are accompanied by a finite rate of spontaneous emission 
scattering, given by T = e J2i m \^i,m\, where e <C 1 is the 
ratio of excited state linewidth to Raman detuning. The 
scaling of these potential errors depends upon the par- 
ticular graph, so we consider two extremes. A uniform 
fully-connected interaction graph can be trivially gener- 
ated with a single spectral component tuned close to the 
COM mode with a detuning \uji -fi\/uJi < logTV/TV 2 . For 
a fixed level of phonon error, the total optical intensity 
should be reduced as logTV/TV, taking into account the 
intensity reduction per ion as the beam is expanded to 
accommodate the linearly expanding chain in space. In 
this case the uniform Ising coupling is expected to scale 
as N\Jij\ oc logTV/TV 2 , and the spontaneous emission rate 
per spin actually decreases with TV. For a sparse interac- 
tion graph, such as a ID (nearest-neighbor) Ising model, 
all modes are involved, and this time for a fixed phonon 
error the total optical intensity can remain fixed, since 
the typical mode splitting falls only as 1/TV, while spon- 
taneous emission per ion is fixed. The calculation in Fig. 
4 shows that the resulting nearest-neighbor interaction 
scales as Jj,i+i oc 1/TV. In either case of fully-connected 
or local Ising model, we thus expect to be able to sup- 
port significant Ising interaction strengths with up to a 
few hundred spins. 

For a general Ising graph, from Eq. 4 we find that 
each pairwise interaction J,- j depends upon a balance of 
TV terms, and errors will accumulate with TV from fluctu- 
ations of relative optical intensities of the various spectral 
components of the beam (which should be stable if the 
spectral components are generated with high quality ra- 
diofrequency sources and modulators as shown in Fig. 
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FIG. 4: Scaling of the nearest-neighbor coupling for the case 
of the production of a ID Ising chain for N — 3 to TV = 33 
ions. The total optical intensity is fixed with J^. „ \^i,n\ = 
1 MHz , and each spectral component is detuned from its 
motional sideband by a fraction f s — 0.03 of smallest mode 
splitting, with a center-of-mass (COM) frequency uj\/2tt — 5 
MHz . We find that the resulting Ising coupling scales roughly 
as l/N. 

3) or their detunings from the motional sidebands. The 
most important source will likely be fluctuations in the 
motional trap frequencies, where we expect the fractional 
error in the Ising coupling to grow as \/N{5u) m /oj m ), so 
that a typical fractional fluctuation in the motional trap- 
ping frequencies of ~ 10~ 3 might be expected to cause 



Ising coupling errors at a level of about 1% for N ~ 100 
ions. 

The scheme presented here can also be applied to more 
general Heisenberg spin models involving other noncom- 
muting spin-spin interactions, such as the XY model or 
the 2D hexagonal Kitaev model relevant to topological 
quantum degrees of freedom [36]. Here, additional Ra- 
man beams that couple to the other axes of motion can 
be exploited [9] . Alternatively, a single direction of mo- 
tion can be used as discussed above, with a stroboscopic 
alternation between Raman laser beams with different 
beat note phases as the various Ising interactions in the 
Hamiltonian are applied sequentially. Here, we employ 
the Trotter expansion of the evolution operator [2, 18] 
and switch the various Ising terms rapidly enough so 
that higher order terms in the expansion can safely be 
neglected. Although this discussion concentrated on a 
linear array of ions in space, these ideas apply in general 
to any stable ion crystal where the motional sidebands 
are resolved and prepared in the Lamb-Dicke limit, and 
should be useful for higher-dimensional trap geometries 
such as trap arrays [23, 24] or Penning traps [11, 37]. 
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